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We present a Heralded Photon Source based only on linear optics and weak coherent states. 

By time-tuning a Hong-Ou-Mandel interferometer fed with frequency-displaced coherent states, the 
output photons can be synchronously heralded following sub-Poisson statistics, which is indicated by 
the second-order correlation function ( g 2 (0) = 0.556). The absence of phase-matching restrictions 
makes the source widely tunable, with 100-nm spectral tunability on the telecom bands. The 
technique presents yield comparable to state-of-the-art spontaneous parametric down-conversion- 
based sources, with high coherence and fiber-optic quantum communication compatibility. 
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I. INTRODUCTION 

Single-photon sources are a fundamental resource in 
quantum optics, with a broad range of applications P- 
0 . Even though deterministic single photon generation 
is highly desirable, it is not yet available. Therefore, 
the photon number statistics of any pseudo-single-photon 
source is always a concern, and the compromise between 
vacuum, single- and multi-photon emission probabilities 
directly impacts the source’s performance. 

A practical and low-cost way to probabilistically pro¬ 
ducing single-photon pulses is to use a faint laser source, 
whose output pulses approach coherent states [5]. Each 
pulse is described as a superposition of Fock states with 
the number of photons following Poisson distribution ac¬ 
cording to the parameter /i, the average number of pho¬ 
tons per pulse. A small value of fi binds the proportion of 
multi-photon to non-vacuum pulses to /i/2 but leads to 
an increase in the probability of emitting vacuum pulses 
to about 1 — fi. The weak coherent states (WCSs) lead to 
a trade-off, as multi-photon and vacuum pulses impose 
limitations to the security and performance of quantum 
key distribution (QKD) links. Safeguarded by strict se¬ 
curity proofs nam and techniques pm however, WCSs 
are widely used in QKD. 

Alternatively to a faint laser source, nonlinear op¬ 
tical processes can be explored to approach a single¬ 
photon source with sub-Poisson photon-number distri¬ 
bution, which is the case of Spontaneous Parametric 
Down Conversion (SPDC) p~TI H3] and Four-Wave Mix¬ 
ing [14j [15]. In such processes, pump photons can be 
converted into pairs of highly-correlated phase-matched 
photons. The detection of one photon can, thus, be 
used to announce the presence of the other photon. This 
so-called heralded single-photon source can be obtained 
from continuous-wave (CW) [T6HT8] or pulsed |19j 20 op¬ 
tical pump even though the inherent low-coherence na- 
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ture of the non-linear processes yield broadband emis¬ 
sion, with ~ 10 ps coherence time m This feature 
renders the synchronization of fiber-optic QKD systems 
critical - specially if it involves the joint measure of two 
single-photon pulses from distinct remote sources (e.g. 
Bell-state measurements [2TU23] ). The coherence time of 
the photons can be enhanced to the ns range by using 
ultra-narrow filtering techniques, for which delicate ther¬ 
mal stabilization is required mm- Nevertheless, the 
narrower the filtering, the less bright is the emission and, 
thus, the throughput of the source. Other rather com¬ 
plex resources, such as double-resonant cavities, could 
also be employed to achieve higher coherence [24]. These 
advanced techniques, however, are frequency-dependent 
providing barely tunable emission. 

In this paper, we propose and experimentally demon¬ 
strate a heralded photon source (HPS) fully based on 
WCSs and completely independent of nonlinear optical 
effects. The Linear-Optic HPS relies on the interfer¬ 
ence of two frequency-displaced WCSs in a Hong-Ou- 
Mandel (HOM) interferometer l 25 j . The interference pat¬ 
tern exhibits anti-bunching peaks with increased correla¬ 
tion between the interferometer’s output spatio-temporal 
modes. The local detection of a photon in one of these 
modes enables the output of the source leading to the 
announcement of a squeezed photon-number pulse. The 
process is based only on linear optics with no phase¬ 
matching restrictions, thus making the proposal freely- 
tunable over a 100-nm wide spectral range on the telecom 
bands without performance loss. 

The paper is divided as follows. Section [IT] presents the 
physical principle driving the proposed source and the 
model describing the photon-number statistics. Section 
mu describes the experimental setup used for generating 
the frequency-displaced WCS and the source implemen¬ 
tation. The Hanbury-Brown and Twiss analysis, a sim¬ 
ulation of a QKD link, and the results are presented in 
Section m Section [V] presents the final remarks and the 
conclusion. 
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II. TWO-PHOTON INTERFERENCE WITH 
WEAK COHERENT STATES 

It has been formerly shown that the incidence of two 
indistinguishable single photons in a beam splitter (BS) 
leads to the photon-bunching effect when the temporal 
modes are matched [25] . A quantum beat pattern is ex¬ 
pected when the photons have different frequencies [[26] 
and the effect can be observed in temporal modes if their 
coherence time is high enough EHI2B]. The optical beat 
note corresponds to the frequency displacement between 
the photons even if independent WCSs are used m • 

Consider a symmetric BS with input spatial modes 
A and B and output spatial modes C and D. The 
model considers that the input states are generated by 
two independent continuous-wave (CW) optical sources 
emitting spatial-mode matched parallel-polarized pho¬ 
tons with identical average number of photons per 
time interval. We decompose the input frequency- 
mismatched WCSs into pairs of Fock states with M 
and N photons, |M, N)a,b- The error associated to 
the simplification of bounding M + N < 3 is esti¬ 
mated from the Poisson cumulative distribution function, 
Em=o ELo H M+N e~ 2 ^/ (M\N\) [5], as being lower 
than 1% for p, < 0.67. 

The case in which any number of photons enter the 
interferometer through one of the input modes, with vac¬ 
uum at the other, |M, 0)a,b or |0, N)a,b, has a straight¬ 
forward solution: a random distribution of photons be¬ 
tween the two output spatial modes. We focus, there¬ 
fore, on the case of two frequency-displaced single pho¬ 
tons entering the HOM interferometer, one in each of 
the input modes. This main result is then used to com¬ 
pose the other non-trivial cases of one- and two-photon 
states entering each input mode 11, 2 )a,b, |2,1 )a,b, and 
|2,2 )a,b- The possible outcomes of the interferometer 
are then combined and weighted by the probability of 
occurrence of the correspondent input states. This in¬ 


formation consists on the source’s output statistics, once 
the probability of a heralding event is considered. 

A. Theoretical Model for the Linear-Optic HPS 

The transfer function for the symmetric BS can be de¬ 
scribed using creation operators 0: 

a A = (j a c + a \ 

V (!) 

°} B = (c^c T J a \^j /v^2 

Combining the description of the wave-packet of each 
input photon £a,b (t) = £a,b (t) e~^ A ’ B to the creation 
operators, we write electric field operators that can be 
associated to the input modes of the BS [28] : 

E A,B (*) = £a,b (t) a A ,B 

E a ,b W = £a,b W ®a,b 

The BS’s output spatio-temporal modes can, therefore, 
be written as a function of the field operators [28] |29] 

E+(t)=[-jE+(t) + E+(t)\/y/2 
E+(t)=[E+(t)-jE+(t) ]/V2 

The case of interest is that of two single photons en¬ 
tering the BS, one at each input mode, which can be 
written in terms of creation operators applied to vacuum 
as |in )a,b = |C 1 )a,b = a \°^ B |0, 0 )a,b- We adopt the no¬ 
tation p rI for the conditional probability of finding 
R and S photons at the output modes C and D , respec¬ 
tively, given that we had M and N photons at the input 
modes A and F>, respectively. 

The probability of finding one photon in C at t = To 
and the other photon in D at t = To -Ft can be computed 
as 


P i,iou t ( T o, r) = (0| a A a B E c (r 0 ) E D (r 0 + r) E £ (r 0 ) E+ (r 0 + r) a ] A a ] B |0) 
resulting in 

-Pl.lout ( T °’ T ) = \ l^( T o + r )£s( r o) - U(to)£b(t 0 + t ) I 2 


Eq. [4] can be solved considering gaussian-shaped wave- 
packets with identical half-width at 1/e, cr, and well- 
defined angular frequency modes uja and uob {u — ^a/ 2+ 
cjb/ 2). The frequency difference A = ujb — oja is fixed 


and the wave-packets can be expressed as 

1 


Cl (*) = 


e -(t-5r/2) 2 /(2a 2 ) e -j(uj-A/2)t 


V 7T<7 

£ b ^ = ^ Q — (t+ST/2) 2 /(Zcr 2 ) c ~j(u+A/ 2 )t 


(5) 


where Sr is the relative delay between the wave-packets. 
Integration from — oo to oc over tq (as we are only 
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interested in the effective time delay r between the wave 
packets) and over Sr (to take into account the CW nature 
of the optical sources) yields [29] 

pl dl = \^-P) ( 6 ) 

r 2 

with /3 = e 2^2 cos (tA). Due to normalization and sym¬ 
metry arguments, the probability of finding two photons 
in one output spatial mode and vacuum in the other is 
found to be 


pldin _ pl,lin 

2,0 ou t 0>2out 


d 1+/3 ) 


(7) 


The solution considering a single photon entering one 
of the input ports of the beam splitter with vacuum en¬ 
tering the other, |l,0)i n and |0,1 )* n , is straightforward: 


pl>0in pl,0in pO,li n pO,li n 1 /O 

^l.Oout _ ^Odout _ *l,0 O ut “ MJ.lout _ i / Z 


(8) 


The same analysis applies to the case in which either 
two (Eq{9| or three photons (EqflO|) enter the same input 
port and vacuum enters the other 


p2,0i n _ p2,0i n _ p0,2i n _ p0,2i n 

^2,0out _ M),2 ou t — ^2,0out MJ^out 

p2,0in _ p0,2 

^MoUt _ "^15 lout 


1/4 

1/2 


(9) 


p3,0in _ p3,0in _ p0,3in _ p0,3i n 

^3,0out _ M),3 0 ut - ^.Oout “ ^0,3out 

p3,0in _ p3,0in _ p0,3in _ p0,3i n 

1 O 1 - JL 1 O - Jl 9 1 - ± -1 9 

^Dout J-^out ^Dout 


1/8 

3/8 


( 10 ) 


The final probabilities for each output state in modes 
C and D are obtained by grouping the conditional prob¬ 
abilities of each output pair, Eqs. weighted by 

the probabilities of occurrence of each pair of input Fock 
states. The latter are given, for WCSs, by the Poisson 
distribution P (M, N\/j,) a B = [5]. The gen¬ 

eral form of the probability of occurrence of the output 
state | R, S)c,d when the BS is fed by frequency-displaced 
WCSs is, thus, given by 


3 3 (R,S,M,Ne N 

Pr,s =EE P R,s Np M,N, l R + S< 3 ( 12 ) 

m= o n=o [ M + N < 3 

B. Statistics of the Linear-Optic HPS 


The output states of the source are conditioned to the 
detection of one - or more than one - photon at the 
spatial mode C, the heralder arm. We thus normalize the 
emission statistics of the source using those cases when 
at least one photon heralds the output mode, which is 
given by 


3 2 


Pt = EE VrPr,s, 

R= 1 S =0 


R,Se N 
R + S < 3 


( 13 ) 


Therefore, the probabilities of the three distinct events 
of heralding vacuum P v , single- P s or multi-photon P m 
pulses at the spatial mode D (disregarding dark counts) 
are given by 


Py = (77lPl,0 + P‘2,0 + V%P 3,0 ) / Pt ( 14 ) 


The cases | 2 , l)i n and 11 , 2 )i n are modeled as a superpo¬ 
sition of a pair of single photons, one in each input port, 
and an independent photon in either one of the input 
ports. 


Pm = VlPl,2/ P T 

(15) 

P s = (t?iPi,i + 772 ^ 2 , 1 ) /Pt 

(16) 


p2 5 lin _ p2,lin _ pl,2i n _ pl,2i n _ _ (1 . 

r 3,0 ou t — M),3out — r 3,0out — ^O^out — g ^ 


0,3out 


3,0 ou t 


0,3 ou t 


■P) 


p2,lin _ p2,lin _ pl,2i n _ pl,2i n _ _ /O _ 0 \ 

^2,l out _ r i,2out _ ^2,lout _ ■ l,2out _ g ^ 


1 ;2 0 ut 


2,1 0 


1 j2 ou t 


where r]i = 1 — (1 — rjcY is the i-photon detection effi- 
(11) ciency of the heralder detector of overall detection effi¬ 
ciency r\c • With a = exp (—2/i), Eqs. ( I4p6 ) can be 
expressed as 


Py = 


P — 

1 m — 


P, = 


24 + n [24 - 12 r) C + (12 - 6 rj c ) /?] + M I 2 [12 - 12 T)c + 4 rfc + (9 - 9r] C + 3 ■rfc) 0] 
24 + Ji (48 - 127/c - GncP) + V 2 [48 - 24 t/c + W c + ( 3 W c - 6 vc) P] 

_ M 2 (12 ~ 3/3) _ 

24 + n (48 - 12 r]c - 6t? C P) + M 2 [48 - 24t? c + Wc + (3 Vc ~ 6 Vc) P] 

M (24 - 12/3) + /j, 2 [24 - 1277c + ( 377 c - 6) /3] 

24 + /x (48 - 12 t?c - 6t?c/3) + M 2 [48 - 24t?c + W c + (3t? 2 - 6t?c) P\ 


(17) 


III. EXPERIMENTAL SETUP 

We report an all-fiber implementation of the Linear- 
Optic HPS source. Two frequency-displaced WCSs cre- 


I 

ate an optical beat note whose anti-bunching peaks 
are explored for heralding Sub-Poisson photon number 
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pulses. The central wavelength of the source is cho¬ 
sen to be 1545.90 nm, with spectral tunability over 100 
nm around the telecom C band. The self-heterodyne 
approach, [29], depicted in Fig. [l] generates the 
two frequency-displaced WCSs up from a tunable laser 
source. The technique avoids stability issues regarding 
both optical power and relative frequency drift. 



FIG. 1. Preparation of frequency-displaced WCSs with a 
self-heterodyne setup. AWG: arbitrary waveform generator; 
FM: frequency modulation; TGR: trigger signal; TLS: tun¬ 
able laser source; VOA: variable optical attenuator; BS: beam 
splitter; OD: optical delay; DG: delay generator; PC: polar¬ 
ization controller; AM: amplitude modulator; CLD: classical- 
level detection unit. 

The attenuated CW signal of an external cavity laser 
diode is frequency modulated with a triangular wave of 
period T = 322.6//S and sent into a balanced Mach- 
Zehnder Interferometer-like setup (MZI). An optical 
switch at one of the MZI’s arms selects the photons with 
constant frequency offset relative to the other arm. The 
optimized duration of the optical switch pulse was exper¬ 
imentally determined to be of 30/rs - 9.3 % of the trian¬ 
gular wave period. This value yields a classical optical 
beat note of 40 MHz with a 3.1 MHz linewidth verified 
by measuring the bright-light version of the optical beat 
note with a p-i-n photodiode and an electrical spectrum 
analyzer at the output of the MZI (Fig. §a). 

The total frequency sweep that characterizes the 
Linear-Optic HPS output linewidth is experimentally de¬ 
termined with a high resolution Optical Spectrum Ana¬ 
lyzer connected to the source’s output (Fig. [2]-b) as 447.5 
MHz. The measured linewidth is square-shaped due to 
the linear sweep of the frequency modulation. The corre¬ 
spondent coherence time of 2.2 ns is comparable to ultra¬ 
narrow filtered SPDC-based HPSs [24] . 

The output of the MZI is converted into an HOM in¬ 
terferometer to build up the Linear-Optic HPS, which 
can be time-tuned so that the local detection at the her- 
alder mode enables the source’s optical output. In order 
to determine the operational point of the source, the in¬ 
terference pattern is measured with a photon statistics 
analyzer - composed of a remote detector triggered by a 
delayed version of the heralding signal - connected to its 
output, as depicted in Fig. [3] All detectors are InGaAs 


o 
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FIG. 2. Measurement of bright versions of the WCSs at the 
classical-level detection unit: (a) beat spectrum and (b) op¬ 
tical spectrum. 


avalanche photodiode single-photon detectors (SPD) op¬ 
erating in gated-Geiger mode with 2.5 ns temporal gate 
width, 15% detection efficiency and 10 fi s deadtime. 



FIG. 3. HOM interferometer with overlapping WCSs and 
photon statistics analyzer. OS: optical switch; C: pulse 
counter. 

Figure [4] depicts the beat pattern of the HOM inter¬ 
ferometer when both the optical switch and the heralded 
pulse’s detection are synchronously scanned with respect 
to the heralding signal. The pattern is characterized by 
the mutual coherence time of the sources and by their rel¬ 
ative frequency mismatch [28] 29]. The theoretical pre¬ 
diction for the normalized coincidence counts |29 

Ccoinc = ^ (2 - cos (tA)) (18) 

is fit to data. 

The operational point for the Linear-Optic HPS is also 
shown in Fig. [4] which is obtained when the optical 
switch is matched and fixed to r = ±7 t/A - the anti¬ 
bunching peak. Under this configuration, only the pho¬ 
tons enabled by the gate opened at the output optical 
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FIG. 4. Interference pattern measured with the photon statis¬ 
tics analyzer for the 40MHz-displaced WCSs with the opti¬ 
cal switch matched to the operational point (blue) and syn¬ 
chronously swept with the detector (black squares). The the¬ 
oretical model is fit to data (red line). 


switch - and thus conditioned to the local detection at 
the heralder detector - are emitted. Even though the 
overall output of our source exhibits a coherence time of 
2.2 ns, the optical switch’s gate width may be a limit¬ 
ing factor for the coherence time. Given the ~ 5 ns gate 
width shown in Fig. [4j however, this does not impact 
on the photon’s coherence. It is worthy noting that, in a 
single-photon interference context, the linewidth is lim¬ 
ited to around 200 MHz according to the Fourier Trans¬ 
form of the output pulse width. 


IV. RESULTS 

A. The Second-order correlation function at zero 
time 

Two distinct values of the relative temporal mode de¬ 
lay, r, are of importance in analyzing the source’s statis¬ 
tics: r greater than the mutual coherence time of the 
WCSs ( r co h ); and r = ±7t/A, the operational point of the 
source which corresponds to the highest anti-bunching 
peak of the interference pattern. In the first case, the in¬ 
terference vanishes since the states are fully distinguish¬ 
able. Any attempt to herald such states will fall into 
the Poisson-like photon-number statistics of a faint laser 
source. 

In the anti-bunching peak case, however, coincident 
counts can reach up to 150% of the distinguishable case 
- see Fig. [4] The normalized heralding probabilities of 
vacuum, multi-photon, and single-photon pulses are pre¬ 
sented below as a function of g and of the heralder SPD’s 
efficiency (rjc) where the value of /3 was approximated to 


— 1 (anti-bunching condition) in Eq. [l7j 

p = 8 + M (4 ~ 2 rjc) + n 2 (l - rjc + rfc/_ 3) 
v 8 + n (16 — 2r] C ) + A* 2 (16 - 6r? c + Vc/ 3 ) 


p - _T_ (20) 

m 8 + (i (16 — 2r]c) + M 2 (16 — 6r]c + Vc/^) 

p = _ 12m + M 2 (10 - 5 rjc) _ 

s 8 + n (16 — 2r] C ) + M 2 (16 - &Vc + Vc/ 3 ) 1 

In Table [IJ we present these normalized heralding prob- 
abilities considering g « 1 and rjc = 1 compared to 
the probability of finding similar pulses in a faint laser 
source. Even though the multi-photon emission proba¬ 
bility is slightly higher, the emission of vacuum is rel¬ 
atively suppressed, and the single-photon pulses occur 
with higher probability. This characterizes the narrow¬ 
ing of the photon-number distribution of the squeezed 
WCS at the output of the source m- 


TABLE I. Normalized Probabilities (g « l,rjc = 1) 


Pulse 

Linear-Optic HPS Faint laser source Ratio 

Vacuum 

1 - \ji + 1 - y + \jj? 1 - \g 

Multi-Photon 

y 2 y 2 i- y 

Single-Phot on 

y - 2/i 2 g — g 2 § - \g 


The second-order correlation function parameter g 2 (r) 
at zero time determines whether the source produces pho¬ 
tons following sub- (< 1), super- (> 1), or Poisson-like 
(= 1) statistics [HI [18]. The experimental values are ob¬ 
tained with a Hanbury-Brown and Twiss (HBT) Ana¬ 
lyzer set at the output of the source. The HBT setup 
consists of a symmetrical BS connected to the heralded 
mode, an SPD connected to each of its output ports trig¬ 
gered by the heralding signal, and a coincidence station 
(CS) (Fig. §. 



FIG. 5. Experimental setup for the HBT Analyzer. 
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The parameter g 2 ( 0) is experimentally defined as 16] 


multi-photon pulse, so we can write 


g 2 ( 0) = Qfg/ ( QfQg ) (22) 


Qf,g 


P s Vf,g 

2 



(23) 


where, Qi?, Qq, and Qfg are, respectively, the proba¬ 
bilities of an event being registered at SPD F, SPD G, 
and in coincidence between them. These values can be 
computed from the probability of an z-photon state be¬ 
ing heralded (the source statistics) and detected in the 
desired detector (the single-photon detection probability 
given by t]f and rjc)- For negligible dark count proba¬ 
bility, a detection event can be triggered by a single- or 


Given our suppositions, coincident detections are trig¬ 
gered by half the multi-photon probability events, result¬ 
ing in 


Qfg = 


VFrjG 

Pm/2 


(24) 


The g 2 (0) parameter is then re-written as a function of 
the photon statistics of the Linear-Optic HPS and the 
detection efficiency values as 


0) 


p s/ 2 + (2 - 1 F /4 - rja/A + VfVg/8) P% l + (2- rj F / 4 - rj G / 4) P s P m 


(25) 


From the modeled single and multi-photon probabil¬ 
ities, we find that g 2 ( 0) is 0.56 for the anti-bunching 
and 1 for the distinguishable cases. Fig. [6] presents the 
predicted values for g 2 (0) and the experimental values. 
The uncertainty in the experimental value found for the 
g (0) parameter comes from the fit of the experimental 
parameters (detection efficiency and average number of 
photons per pulse) to the measured gain (obtained with 
the setup shown in Fig. [5| and from the truncation error 
on the model due to the three-photons approach. 



FIG. 6. Theoretical values and experimental results for g 2 (0) 
at the operational point. Inset shows the g 2 (r) function. 

The inset exhibits the predicted values for g 2 (r) con¬ 
sidering g = 0.1 and the detectors efficiency t]h,f,g = 
0.15. The agreement between experimental and theoreti¬ 
cal results, shown in the main figure, validates the model 
and indicates the sub-poissonian character of the source 
due to g 2 (0) < 1. 


B. Simulation of a quantum key distribution link 

The performance of a QKD system is directly de¬ 
pendent on the photon statistics of the optical source. 
The secret key generation probability and the maxi¬ 
mum achievable link distance are analyzed below for 
our Linear-Optic HPS, for a faint laser source, and for 
an SPDC-based HPS, following both the GLLP security 
analysis [ 6 ] and the decoy states approach [ 8 HT 0 ]. The 
parameters Po, Pi and P 2 are the probability of Alice 
sending vacuum, single- or multi-photon pulses for each 
kind of source considered. 

The yield is the conditional probability of detection at 
Bob, given that Alice has sent an z-photon pulse. This 
is given by Yi « Pd + r\u where Pd is the detector dark 
count probability and rji = 1 — (1 — 77 ) L The overall 
transmittance g is composed by the product of the ef¬ 
ficiency of Bob’s detection devices ( r\Bob ) and the link 
loss ( 10 -aL / 10 ) - where a is the attenuation coefficient 
(dB/km) of the fiber with length L (km) at the operation 
wavelength. The gain of the z-photon state is given by 
Qi = PiYi , and the overall gain is obtained by summing 
over the contribution of all states Q p = Qi- 

The error probability of the z-photon state is given by 
Ci a* (eoLo + e opt r]i ) /Y$, where eo=l /2 is the probability 
of a dark count to occur at the wrong detector and e opt 
corresponds to the optical misalignment of Bob’s appa¬ 
ratus. The overall error is given by E p = e iQi- 

The secret key generation probability is computed as 

R>q{Qi[l- H 2 ( ei )] - Q p H 2 (E p ) f (P M )} (26) 

where H 2 (x) is the Shannon binary entropy, f{pc) is the 
inefficiency of the error correction (1.16 here), and q=l /2 
appears due to the bases matching probability. 
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Considering the GLLP security analysis [6], the gain 
of single-photon pulses is lower bounded by 


Qi = Q l *-'52Pi™Q li -P m {n) (27) 

i =2 

This assumption is the pessimistic one where all multi¬ 
photon pulses sent by Alice are eavesdropped. In this 
analysis, the single-photon error is given by 

ei = E ll QJQ 1 (28) 

When using the decoy states method, the values of 
Qi and e\ can be estimated without the pessimistic as¬ 
sumption. This is possible due to the random choice of 
intensities sent by Alice. The i-photon yield does not de¬ 
pends on Alice’s choice, so Eve cannot predict nor fake 
the output of Alice’s source. Furthermore, due to these 
tighter bounds, Alice can usually use higher intensity val¬ 
ues without jeopardizing the system security. 

The simulation parameters were extracted from 
[S]: ck=0.21 dB/km; rj B ob=0M5; P d = 0.85xl(T 6 ; 

e opt =0.033. The value of g was optimized at each link 
distance for the faint laser source and Linear-Optic HPS. 
The efficiency of the her alder detector of the Linear-Optic 
HPS is 15%. The probabilities for the SPDC-based HPS 
were obtained from [19]: Pi =0.42 and g 2 (0)=0.018, with 
P 2 = g 2 (0) Pf/2 m- 

The results for the secret key generation probability 
obtained as a function of the link length for different 
sources are shown in Fig. [7] 



FIG. 7. Simulation of the secret key generation probability for 
a BB84-based QKD system considering: GLLP analysis with 
(a) faint laser source, (b) SPDC-based HPS, and (c) Linear- 
Optic HPS; decoy states method with (d) faint laser source, 
(e) SPDC-based HPS, and (f) Linear-Optic HPS; and (g) a 
true single-photon source. 


C. Discussion 

The reported Linear-Optic HPS, when conditioned to 
the temporal selection of pulses at the operational point, 
exhibits improved statistics when compared to that of 
coherent states. The output states may be regarded 
as squeezed photon-number distribution weak coherent 
states since the second-order correlation function at zero 
time is less than one 0130]. 

The present QKD link simulation reveals that, under 
GLLP analysis, the proposed source increases the max¬ 
imum achievable distance of a QKD link with a factor 
of 1.15 when compared to both the faint laser source 
and SPDC-based HPS. Considering short links, the se¬ 
cret key generation probability of the source is (at least 
1.86 times) greater than for the faint laser source, while 
largely overcome by the SPDC-based HPS (with a factor 
i6 ). ‘ 

Considering the decoy states approach, we found that 
the proposed source outperforms the faint laser source 
in secret key generation probability and in distance. Its 
performance is also highly competitive with state-of-the- 
art SPDC-based HPS concerning both parameters. The 
maximum achievable distance is slightly shorter (4 km, 
corresponding to a length reduction of 2.5%) and the 
secret key generation probability presents a small penalty 
on the simulation. Parameters taken from were also 
tested with qualitatively similar results not shown in the 
figure. 

V. CONCLUSION 

We have presented for the first time, to our knowl¬ 
edge, a heralded photon source relying on linear optics 
and weak coherent states. The interference between the 
frequency-displaced weak coherent states in an Hong-Ou- 
Mandel interferometer results in a structured interference 
pattern that exhibits anti-bunching peaks. The system 
is time-tuned so that only photon pulses that match the 
highest of those peaks are heralded. This enhances the 
correlation between the locally-detected spatio-temporal 
mode and the heralded mode. Sub-Poissonian photon 
statistics can be synchronously achieved. The photon 
statistics of our linear-optic heralded photon source was 
modeled and an Hanbury-Brown and Twiss analysis as¬ 
sessed its sub-Poisson character, with a second-order cor¬ 
relation parameter down to 0.556. 

The performance of the Linear-Optic HPS was evalu¬ 
ated in a simulation of a BB84-like QKD section. The 
source was compared to a faint laser source and to 
an SPDC-based heralded photon-source under both the 
GLLP security analysis and the decoy states approach. 
The simulation results reveal that the proposed source 
outperforms the faint laser source and is highly compet¬ 
itive with the SPDC-based HPS. 

The present 2 kHz heralding rate of the self-heterodyne 
scheme can be improved by employing two continuous 









wave lasers with fixed displaced frequencies. We estimate 
that the heralding rate can be as high as 1.5 MHz, a rate 
comparable to state-of-the-art SPDC HPSs [20] . This 
value is achieved when we assume a 100 MHz triggering 
rate for the local detector, with a detection efficiency 
of 15%, and CW independent laser sources emitting 0.1 
photons per time interval. 

The absence of phase-matching restrictions, as linear 
optics is used, renders the proposed source freely-tunable 
over a 100-nm wide spectral range in the telecom bands 
without performance lost. The physical principle could 
also be explored to assemble a heralded photon source in 
the visible spectral range, matching free-space quantum 


communications. 

In conclusion, the all-fiber telecom-compatible design, 
its sub-poissonian character, and the rather simple setup 
make the Linear-Optic Heralded Photon Source a proper 
choice for long-distance quantum communication, broad¬ 
ening the plethora of technological resources in the field. 
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